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Abstract. We study the square-lattice XY model in the presence of random phase 
shifts. We consider two different disorder distributions with zero average shift 
and investigate the low-temperature quasi-long-range order phase which occurs for 
sufficiently low disorder. By means of Monte Carlo simulations we determine several 
universal quantities which are then compared with the analytic predictions of the 
random spin-wave theory. We observe a very good agreement which indicates that the 
universal long-distance behavior in the whole low-disorder low-temperature phase is 
fully described by the random spin- wave theory. 
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1. Introduction 

The two-dimensional XY model with random phase shifts (RPXY) describes the 
thermodynamic behavior of several disordered systems, such as magnetic systems 
with random Dzyaloshinskii-Moriya interactions Josephson junction arrays with 
geometrical disorder [2, |3], crystal systems on disordered substrates pfl|, and vortex 
glasses [5] . See [6j [7] for recent reviews. The RPXY model is defined by the Hamiltonian 



where ip x = e l9x , U xy = e lAxy , and the sum runs over the bonds (xy) of a square 
lattice. The phases A xy are uncorrelated quenched random variables with zero average. 
A Gaussian distribution, 



has been considered in most of the studies of the RPXY. We denote the RPXY with 
Gaussian-distributed phases by GRPXY. The pure XY model is recovered in the limit 
o~ — > 0, while the so-called gauge glass model [5] with uniformly distributed phase shifts 
is obtained in the limit a — > oo. 

The pure XY model shows a high-temperature paramagnetic phase and a low- 
temperature phase characterized by quasi-long-range order (QLRO) controlled by a line 
of Gaussian fixed points. In the latter phase, the spin-spin correlation function (ip x ipy) 
decays as l/\x — y\ v ^ T ' for large \x — y\, with r\ depending on T; i] ~ T for small values 
of T. The two phases are separated by a Kosterlitz-Thouless (KT) transition [9] at [10] 
Pxy = ^-/Txy = 1.1199(1). For r = T/T X y — 1 ^ + the correlation length diverges 
exponentially as ln£ ~ t -1 / 2 and the magnetic susceptibility behaves as x ~ £ 7//4 > 
corresponding to r\ = 1/4. 

In this paper we shall discuss the low-temperature behavior of RPXY models for 
small disorder. It has already been investigated in several works, most of them focussing 



on the GRPXY; see, e.g., [H El El HU [12l El [13l [HI UHl El [161 HEl UHl EQl EU E21 E31 

ll[2i[2g[iH[27i[iii[iH[iB 



The expected T-a phase diagram, which is sketched in Fig. [H presents two phases at 
finite temperature: a paramagnetic and a QLRO phase. The paramagnetic phase is 
separated from the QLRO phase by a transition line, which starts from the pure XY 
point (denoted by P in Fig. CQ at (a = 0, T = Txy ~ 0.893) and ends at a zero- 
temperature disorder- induced transition denoted by D at (cr ,T = 0). Note that QLRO 
is observed only up to a maximum value o~m of the disorder parameter, which is related 
to the point M = (a M ,T M ), where the tangent to the transition line is parallel to the 
T axis. The transition line from M to D is believed to run (almost) parallel to the T 
axis, with a < a M ; see, e.g., (6J. [j] 

| The first renormalization-group (RG) analyses based on a Coulomb-gas representation of the 
models [I] predicted reentrant transitions for any value of a < n/8, so that <r n = 0. It was then 



H = RevplUijipy cos(9 x - 6 y - A xy ), 



(1) 





(2) 
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Figure 1. Generic phase diagram of two-dimensional RPXY models in the 
temperature-disorder plane. The full line corresponds to the transition line. The 
dashed line is the N line T = a and refers only to the CRPXY model. 



The QLRO phase of the pure XY model is expected to survive for sufficiently 
small values of o (see, however, [38] for a recent critical discussion of this scenario). It 
disappears for large disorder, for [1] a ^ <tju ~ tt/8. The RPXY model for large disorder, 
and in particular the gauge-glass model (<r — > oo), has been much investigated [5] [8] 

tiiiiiaiiiaiisiisEiEa no 

long-range glassy order can exist at finite temperature [22] E3] . Some numerical works 
support a zero-temperature transition; see, e.g., [4"5 ] |4"3] 132"] . 

Beside the GRPXY, we also consider a RPXY model with a slightly different 
distribution, given by 

P(A xy ) oc exp ( . (3) 



V ° 

We denote the RPXY model with the distribution © by CRPXY. The CRPXY is 
interesting because it allows to obtain exact results along the so-called Nishimori (N) 
line 



T = a, (4) 

exploiting gauge invariance [IE] ESI US]- For example, the energy density is exactly 
known along the N line: E = —I\{0)/Iq{0), where (3 = 1/T and I n (x) are the modified 

clarified that this was an artefact of the approximations. Indeed, in agrement with the experimental 
findings pi], numerical works [TT1 [12 H31 IHD] and more careful RG analyses [HI EH EH EH EH [35] 
showed the absence of reentrant transitions for sufficiently small values of a, that is predicted <7o > 0. 
For the CRPXY, [19] suggested that the transition line is parallel to the T axis below the Nishimori 
point M. It is not clear whether this conjecture is correct. Indeed, the analogous conjecture fails in 
the case of the 2D ± J Ising model; see [47] and references therein. 
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Bessel functions. Another important feature of the CRPXY is that along the N line the 
spin-spin and overlap correlation functions are equal, 

where the angular and square brackets indicate the thermal average and the quenched 
average over disorder, respectively. The phase diagram is expected to be analogous to 
that of GRPXY model. As already noted in [4j5], the N line must play an important 
role in the phase diagram, because it is expected to mark the crossover between the 
region dominated by magnetic correlations and that dominated by randomness. In [19] 
it was proven that the critical value % of a along the N line is an upper bound for 
the values of a where magnetic QLRO can exist (note that this does not exclude the 
existence of a glassy QLRO for o > o~m). Therefore, we can identify it as the point 
M where da c /dT = in Fig. [TJ The critical value a at T = must satisfy cr < a M , 
leaving open the possibility of reentrant transitions. 

In this paper we focus on the QLRO phase of RPXY models, and investigate its 
nature. In particular, we check the random-spin- wave scenario pQ, in which the long- 
distance behavior is essentially identical to that in the model obtained by replacing 

cob{6 x -9 v - A xy ) — ► 1 - 1 -{Q X - 9 y - A xy f (6) 

in the Hamiltonian ([1]). We numerically check this scenario by perfoming Monte Carlo 
(MC) simulations of the GRPXY and CRPXY models and provide conclusive evidence 
of the existence of a low-T phase with QLRO determined by the random spin-wave 
theory. Some numerical evidence of QLRO in the GRPXY model was already presented 
in [30]. 

The paper is organized as follows. In Sec. [2] we summarize the main predictions of 
the spin-wave theory, which are then compared with the numerical data in Sec. El Our 
conclusions are presented in Sec. HI The definitions of the quantities we consider are 
reported in Appendix A In Appendix B we report the spin-wave calculation of r\ s in 
the CRPXY model. 



2. The random spin- wave theory 

In the spin-wave limit the partition function is given by 

Z{{A}) = I [d<P]e- H ^ T , H sw = 1 - h ~ ^') 2 ' ( ? ) 

J (v) 
where the link- variables Aij are uncorrelated quenched random variables with Gaussian 
probability P(Aij). For the GRPXY the spin and the overlap correlation functions are 
given by [11139"] 

G s (x -y) = [(e^-W)] = exp[(T + a)G(x - y)], (8) 
G (x -y) = [|< e ^-«)| 2 ] = exp[2TG(x - y)}, (9) 
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where G(x — y) is the (infrared-regularized) two-point function of the Gaussian model 
without disorder 

For the CRPXY one should take into account the nontrivial dependence of P(Aij) on Ay. 
For the overlap correlation function one still obtains for any probability distribution 
G (x — y) does not depend on randomness in the spin-wave approximation. For the spin 
correlation function the a dependence at T = is more complex. For o — > we obtain 

G s (x -y) = exp[(T + a + a 2 /2)G{x - y)], (11) 

disregarding terms of order a 3 . The derivation is reported in|Appendix B 



The above-reported results allow us to evaluate the exponents r] s and r] which are 
related to the large-distance behavior of the spin and overlap correlation functions: 

[ (GRPXY) 

T + ff + a/2 (CRPXY) 
^ 2ix 

G (r)~r-*, r]o = -- (13) 

7T 

We can thus write for both models, at this level of approximation, 

G s (x -y) = exp[2nr] s G(x - y)\ G Q {x-y) = exp[2nr] G(x - y)]. (14) 

Note that the disorder dependence is completely included in the exponents n s and r\ . 
Using these expressions we can campute the universal functions R s {rj s ) and R (rj ), which 
express the ratios R s = and R Q = £ /L in terms of the corresponding exponents 
i] s and rj . It is clear that R s (r] s ) and R (rj ) are identical [R s (x) = -R D (^)] and disorder 
independent, hence they coincide with those relevant for the pure XY model. These 
functions are shown in Fig. [2j Below we show numerically that these predictions are 
satisfied by our numerical data. This provides a clear evidence for the spin-wave nature 
of the QLRO in the RPXY models. A similar strategy was applied in [50] to clarify the 
nature of the low-temperature phase of fully frustrated XY models. 



3. Monte Carlo results 



In this section we numerically investigate the nature of the QLRO phase that occurs 
for sufficiently small disorder. In particular, we want to provide a stringent check 
of the random spin-wave scenario. The quantities which we compure are defined in 
Appendix A 



3.1. Numerical details 

We performed MC simulations of the GRPXY and CRPXY models, considering square 
lattices of linear size L with periodic boundary conditions. We set in both cases 
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Figure 2. R s = £ S /L versus T] s (above) and R = £ /L versus r\ (below). The full 
lines are the results of the random spin-wave theory, obtained by using the expressions 
(fT4l) . Th e MC data are obtained for the CRPXY and GRPXY models with a = 0.1521. 

a = 0.1521, which is well below the maximum value <ta/ [om ~ 0.30 in the CPRXY 
model [51] and o M ~ tt/8 in the GRPXY model) and considered several values of T 
below the critical temperature T c (a), which marks the end of the QLRO phase. MC 
simulations in the high-temperature phase [51] indicate T c = 0.771(2) for the GRPXY 
and T c = 0.763(1) for the CRPXY. 

In the simulations we used a mixture of standard Metropolis and overrelaxed 
microcanonical updates: A single MC step consisted of five microcanonical sweeps 
followed by one Metropolis sweep. In the MC simulations of the CRPXY model we 
also used the parallel-tempering method [521 [53]. It allowed us to perform efficient 
simulations in the region T < a. In the parallel-tempering simulations we considered 



□ CRPXY 
o GRPXY 
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Figure 3. MC estimates of \iiXs/L 2 versus \nL for the CRPXY model at a = 0.1521 
and two values of (3 = 1/T, corresponding to T w 0.139 and T sa 0.714. 



At systems at the same value of a and at Nt different inverse temperatures (3 m - m = 
■ ■ • , /5 m ax- The largest value /3 max corresponded to the minimum temperature value we 
were interested in. The value (3 min was chosen in the paramagnetic phase and was such 
that thermalization at /3 — f3 m \ n was sufficiently fast. The intermediate values j3i were 
chosen such that the acceptance probability of the temperature exchange was at least 
5%. Morever, we always included the value /? = 1/cr, which lies along the N line. This 
provided a check of the numerical programs, since the MC results could be compared 
with the known exact results [H]. The overlap correlation functions and corresponding 
Xo arid £ were obtained by simulating two independent replicas for each disorder sample. 

In the case of the GRPXY model we performed standard MC simulations at 
T = 2/3, 1/2, 2/5, for lattice sizes 10 < L < 40. Typically, we considered 50000 disorder 
realizations and performed O(10 5 ) MC steps for each of them. In the case of the CRPXY 
model the parallel-tempering method allowed us to investigate the temperature range 
T < T c down to T k 0.139, which is below the N line, i.e. satisfies T < a = 0.1521. 
We performed simulations for lattice sizes 10 ^ L < 30. Typically, we considered 25000 
disorder realizations. For /3 = 1/T = 7.2 we also performed standard MC runs up to 
L = 85. 



3.2. QLRO in RPXY models 

We estimate the exponents r] s and r] by studying the finite-size behavior of the 
susceptibilities Xs and Xo defined in Appendix A Indeed, for L — ► oo they behave 
as 

Xs,o - L 2 -^° . 



(15) 
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Figure 4. The MC estimates of rj s and r/ a vs T for the GRPXY and CRPXY models 
at (7 — 0.1521. The lines shows the spin-wave approximations (fT2|) . and (fl3l) . The two 
lines that give rj s for the GRPXY and CRPXY models cannot be easily distinguished 
on the scale of the plot. The dotted vertical line corresponds to T = a = 0.1521; in 
the CRPXY model this point belongs to the N line. 



Estimates of Xs are shown in Fig. [3j On a logarithmic scale, the data fall on a straight 
line, indicating that the asymptotic behavior ( fl5l) already holds for the values of L we 
consider. In Fig. H] we show the estimates of i] s and n for the GRPXY and CRPXY 
model. For T < 0.2 they agree with the spin-wave approximations ( Tl2|) and ( fl3l) . 
Morever, they appear to be mostly independent of the model, in agreement with the 
random spin- wave predictions (note that cr 2 /2 = 0.0116 for a = 0.1521). For a more 
quantitative check, in Fig. [5] we plot the difference 2n s — r] vs T, and compare it with 
the low-order approximations 



2Vs -Vo = ~ (GRPXY), (16) 

7T 

2vs -Vo= ° + (CRPXY). (17) 

7T 

The agreement is very good. Morever, the above-reported relations appear to hold up 
to temperatures close to the KT transition T c (for a = 0.1521 we have T c = 0.771(2) for 
the GRPXY model and T c = 0.763(1) for the CRPXY model), suggesting that they may 
also hold at the KT transition. Since MC simulations in the high-temperature phase [JT] 
show a clear evidence that r] s = 1/4, this may suggest that at the KT transition 

Vo « I - - (GRPXY), (18) 

1 a + a 2 /2 tn 
Vo^t; (CRPXY), (19) 

Z 7T 

The most important check of the spin-wave nature of the QLRO is provided by the 
MC data shown in Figs. [H where we plot R s vs r] s and R Q vs r\ : they agree with high 
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Figure 5. We plot the difference 2i] s — r\ vs T at <r = 0.1521, and compare it 
with the low-order spin-wave approximations for the GRPXY and CRPXY models, 
cf. (|16p and (|17[). respectively dotted and dashed lines. The vertical dotted and dashed 
lines indicate the critical temperatures of the two models, i.e. T c = 0.771(2) and 
T c = 0.763(1) for the GRPXY and CRPXY, respectively. 



accuracy with the curves R s {rj s ) and R (t]o) computed in the spin- wave limit. We believe 
that these results provide a conclusive evidence that the QLRO phase is determined by 
random spin-wave theory. 

Finally, we report some results for the gauge-invariant correlation function (1A.4[) . 
see Appendix A In the spin-wave limit one finds jl9] 

[( e *W*-Ay— •Vv-W)] = exp[(T - a)G(x - y) - \x - y\a/2], (20) 

which predicts that gauge-invariant spin-spin correlation functions are not critical. 
For instance, in the large-L limit the correlation length £g defined from the large- 
distance exponential decay of the gauge-invariant correlation function is finite and given 
by £^ gap ) = 2/cr, independently of T. These predictions are confirmed by our MC 
simulations. We compute the second-moment correlation function £ g defined in (1A.5jl . 
The results are reported in Fig. [6j It is evident that £ 9 is finite in the large-L limit, 
satisfies £ 9 ^ 2/er ~ 13, and is independent of T. Again, this result shows that the 
critical behavior is correctly described by the spin-wave theory. 



4. Conclusions 



In this paper we have studied the low-temperature low-disorder phase in RPXY models. 
We have considered two different disorder distributions and for each of them we have 
computed numerically the exponents rj s and T] , and the correlation lengths £ s and £ G . 
These results have been compared with the predictions of the random spin-wave theory. 
Our main results are the following. 
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Figure 6. The correlation lengths £ s , defined from the standard spin-spin correlation 
function, and £ g defined from a gauge invariant spin-spin correlation function, versus 
L. Results for the GRPXY model with a = 0.1521. 

1) The ratios C, s /L and £, /L, when expressed in terms of the corresponding exponents 
t] s and 7] , are in perfect agreement with the analytic predictions. This indicates 
that the expressions fHM hold quite precisely in the whole low-temperature QLRO 
phase . 

2) The expressions ffl2l and ffl3l) hold only for very low values of T. However, the 
difference 2r] s — r] is apparently well described by spin- wave theory up to the critical 
transition which marks the end of the paramagnetic phase. 

3) In agreement with the random spin-wave theory, the gauge-invariant spin 
correlation function ( 1A.4I) is not critical. 

Finally, note that our calculations refer to probability distributions for which [A] = 0. 
Very little changes if we consider a nonzero average; for instance, one might use the 
distribution 

P(A) oc exp [-(A - a) 2 /2a] . (21) 

In this case we have [A] = a. The new model can be mapped into the original one by 
performing the gauge transformation 

Vkxa) = e~ ia(Xl+X2) ^ 4* = - a- ( 22 ) 

Hence this model has the same phase diagram as the original one. The trasformation 
( f22l) leaves the overlap correlation functions unchanged, since they are gauge invariant. 
The behavior of the spin-spin magnetic correlation functions is more subtle, lib = (a, a), 
in Fourier space we have 

G s {q-a) = G s (q + b;a = 0), (23) 
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where G s (q; a) is the Fourier transform of the spin magnetic correlation function in the 
theory with a nonvanishing average a. In the standard theory (a = 0) the critical modes 
are those with q = 0, while for q ^ the behavior is not critical. This implies that 
the critical modes in the theory with a ^ are those associated with a nonvanishing 
momentum q = —b. Hence, in this theory the magnetic susceptibility, which corresponds 
to q — 0, is not critical. 

As a final comment, we note that the distribution functions of the phase shifts 
are not gauge invariant: Phase shifts that differ only by a gauge transformation have 
different probabilities. Another interesting issue is whether the results for the QLRO 
phase reported here also apply to gauge-invariant distributions. 



Appendix A. Notations 

In terms of complex site variables ipi = e %e \ the RPXY Hamiltonian takes the form 

(ij) 

where = e tAiJ . 

We consider several two-point correlations functions: the magnetic spin-spin 
correlation functional 

G s (x - y) = Re WM\ = WM\, (A.2) 
and the overlap correlation function 

G {x-y) = [\{rM\\ (A.3) 

which can be written as G Q (x — y) — [(q x q y )}, where q x = ipx ipx and the 
upperscripts refer to two independent replicas with the same disorder. The angular and 
square brackets indicate the thermal average and the quenched average over disorder, 
respectively. We also consider a gauge-invariant spin-spin correlation function 

G g (x-y) = [Re(r x U[T X]y }^ y )}, (A.4) 

where T x . y is a path that connects sites x and y and ^[r^] is a product of phases 
associated with the links that belong to T x . y . The paths connecting the points x and y 
are chosen along the lattice axes, choosing the shortest path (see [M] for details). 

We define the corresponding susceptibilities: the magnetic susceptibility Xs = 
J2 X G s (x), the overlap susceptibility Xo = J2 x Go(x), and Xg = J2 x Gg( x )- We also 
define the corresponding second-moment correlation lengths 

c2 = G#(0) - G # (q min ) 

t /min l - T #vymmJ 

where g min = (2tt/L, 0), q = 2 sing/2, and G#((?) is the Fourier transform of G#(x), and 
# indicates s, o, g. 



The last equality in (|A.2|) can be proved by using the symmetry ip w — > ip* and U xy — > U, 



■j:y 
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In this Appendix we wish to derive (1111) . We follow closely [lj. We first consider the 
spin-spin correlation function G s (r,A a ) = (exp [i(4>(0) — <fi(r))]) for fixed values of the 
random phases A a . As in [TJ we rewrite it as 

(exp [i(flQ) - 0(r)) - Pjd 2 rA ■ V0]) o 
s[r) (exp[-f3jd*rA-V<p}) ' 1 ' 

where (-)o indicates the average with Hamiltonian 

H = | J d 2 r(V0) 2 . (B.2) 

Repeating the steps discussed in [1] we end up with 

G s (r, A) = exp (rG(r) + \ j d 2 s^ A a (s)M a (s, r) j , (B.3) 

where G(r) is the Gaussian propagator ([TO]) and 

M a (s,r) = 2 J -0-^^-1)^. (B.4) 

Note that M a (s,r) is imaginary, M Q (s, r)* = — M a (s, r), so that 

|G s (r,A)| 2 = e 2TG « . (B.5) 

Thus, irrespective of the phase distribution, the overlap correlation function does not 
depend on o. To compute the spin correlation function we must average G s (r,A) over 
the distribution of the phases A a . We consider the general distribution 

P(A) oc exp (-Q^p-) , (B.6) 

which satisfies Q(z) = z for 2: — > and is such that, for o — > the distribution is peaked 
around A = 0. Thus, to compute the expansion of G s (r) for small cr, we can expand 
Q(A 2 ) in powers of A 2 . We assume 

Q(z) = z + az 2 + 0(z 3 ), (B.7) 

where a is a distribution-dependent coefficient. For the distribution (jHJ) we have 
a = —1/12. To compute the correction of order a 2 to r] s we rewrite ([-}a indicates 
the average over A) 



S = 



oc 



exp ^ J d 2 s^A a (s)M a (s,i 
y [DA] exp (~ j d 2 s[A 2 + a(A 2 ) 2 ] + ^ j d 2 sJ2^ a M a j .(B.8) 



We introduce a new field -B a defined by 



A a = V^B a + ~Af a (B.9) 
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and perform the integral over B. Disregarding terms of order a 3 we obtain 



S = exp |(l-6a<r) J d 2 sM(s,r) 2 . (B.10) 

Since 

r d 2 sM(s,r) 2 = 8G(r), (B.ll) 
we obtain finally 

G s (r) = e (^+-6«- 2 )G(r)_ (B.12) 
If we set a = —1/12, we obtain result ffTTT) . 
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